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Abstract. We present one-dimensional KKR, method with the aim to elucidate 
its linear features, particularly important in optimizing the numerical algorithms 
in energy bands computations. The conventional KKR equations based on the 
multiple scattering theory as well as novel forms of the secular matrix with 
nearly linear energy dependency of the eigenvalues are presented. The quasi- 
linear behaviour of these eigenvalue functions appears after (i) re-normalizing the 
wave functions in such a way that 'irregular' solutions vanish on the boundary 
of the 'muffin-tin' segments and (ii) integrating the full Green function over the 
whole Wigner-Seitz cell. In addition, using the aforementioned approach we derive 
one-dimensional analog of the generalized Lloyd formula. 

The novel KKR approach illustrated in one-dimension can be almost directly 
applied to the higher dimensional cases. This should open prospects for the 
accurate KKR band structure computations of very complex materials. 



1. Introduction 

In the Korringa-Kohn-Rostoker (KKR) method the main problem in performing 
electronic structure calculations is due to strongly nonlinear dependence of the 
KKR- matrix eigenvalues on energy for a fixed fc-point. It is very difficult to find 
properly all the zeros of eigenvalue functions, when there is a large number of atoms 
in the unit cell. Hence the fast and accurate extraction of all zeros of eigenvalues still 
remains a challenge in computation of complex materials. 

Butler pQ and Schwitalla et al. [2] formulated the KKR formalism for one 
spatial dimension, which is very similar to its three-dimensional counterpart. The 
one-dimensional case is computationally much simpler and allows to formulate the 
substantial results in compact, analytical forms. Therefore, it is very attractive to 
investigate future of the KKR method along this line. 

In this paper we use the Green function approach for solving the Schrodinger 
equation. As already known, the one-dimensional (ID) KKR method gives exact 
solutions. This is due to the facts that (i) the unit cell can be entirely filled by 
the non-overlapping (touching each other) 'muffin-tin' segments and (ii) the multipolc 
expansion contains only two 'spherical harmonics' (I — and I = 1). Furthermore, the 
KKR structure constants in the ID case can be obtained analytically, which markedly 
facilitates numerical calculations. 

Our paper is organized as follows. In Sec. 2 we recall after Butler pQ and 
Schwitalla et al. [2] the standard 1D-KKR formalism. In this section we also discuss 
few difficulties in determining the band structure in multi-atom systems if using the 
standard approach. Sec. 3 presents the novel form of the KKR-matrix in terms of 
logarithmic derivatives as well as the advantages of this approach with respect to the 
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standard KKR theory. Next, the generalized Lloyd formula 6 is derived , which helps 
in computation of the total number of states and placing properly the Fermi energy at 
early stage of calculations in real materials. In Sec. 5 we show illustrative results for 
an arbitrary chosen potential, which are followed by summary and conclusions (Sec. 
6). 

2. Standard KKR formalism in one dimension 

We look for a solution of the ID Schrodinger equation with Hamiltonian (atomic units 
are used, energy unit is 1 Ry): 

H = ~ + V(x). (1) 
ax 

The potential in Q is assumed to be translationally invariant with the period of lattice 
constant a, so it can be written as 

V(x) = v{x-na). (2) 

n— — oo 

We assume v(x) in J5J to vanish for \x\ greater than some radius S, and do not 
overlap each other. This assumption has the form analogous to the non-overlapping 
'muffin-tin' potential in three-dimensional (3D) KKR theory. We shall expand the 
potential v(x) in @ and the wave functions around each center into symmetric and 
antisymmetric functions, analogous to the 3D expansion into the spherical harmonics. 
First, let us define ID spherical coordinates centered at na-point, 



x e 



a a 
— , na H — 
2' 2 



(3) 



with discrete directional coordinate, x =sgn(x), and radius r £ [0, |] . Then, with the 
definition of ID analog of the spherical harmonics 

Y (x) = 1/V2 Y x {x) = x/y/2, (4) 

we decompose v(x) into spherical and aspherical parts: 

1 x 

v (x) = -^[v(r)+v{-r)} + - [v(r)-v(-r)] = v (r)Y (x)+v 1 (r)Y 1 (i).(5) 

Similarly, we may write ID multipole series for the wave function 

i)(x) = Mr)Yo(x) + Mr)Yi(x) (6) 

with the radial functions of the free Schrodinger equation given by the following regular 
jo(VSr) = cos(\/£V) ji(VEr) = sin(v A Er) (7) 

and 'irregular' 

no{s/~Er) = s'm(VEr) ni(V~Er) — — cos(V~Er) (8) 

special solutions. Also analogs of the 3D spherical Hankel functions, usually defined 
as hi(x) — ji(x) + ini(x), exist: 

h (V~Er) = exp(i\f~Er) hi(VEr) = —i exp(iV~Er) . (9) 
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The square root v E is taken on the complex energy plane with the branch cut along 
the positive real axis. To proceed further we follow the steps in the multiple scattering 
formalism [HI Ej , leading to the following form of the Bloch- Fourier Green's function 

< x'\G(E,k)\x >= Ji{E, X> )Zi(E )X< ) 

(10) 

+ Ei'jZv(E,x') [t- l {E)-B{E,k)] n Z x {E,x) . 

In (|10(l Zi(E, x) and Ji(E, x) are the regular and irregular solutions of the Schrodinger 
equation with the radial parts given in the asymptotic region (r > S) by 

' ^ - (ID 



[Z(E,r)]i>i =jv{VEr)[r 1 {E)] vl - - 7 =hi{VEr)5 V i 

V E 



[J(E,r)] in =ji>(VEr)6m 



(12) 



B{E, k) is the matrix of KKR structure constants with explicit form generalized to 
multi-atom case as given in Appendix. The matrix t(E) is one-scatterer i-matrix and 
may be found from the wave functions 



iPl(E,x) = J2Yi'(Z)bP(E, 



r in 



with radial solutions satisfying the following set of differential equations 
d 2 



dr 2 



ip(E, r) + [E- v(r)]tp(E, r) = 0, 



with the potential v(r) in the matrix form 

v g (r), v a (r) 



v(r) 



v a {r), v s (r) 



1 

V2 



v (r), vi(r) 
Vx[r), v (r) 



(13) 



(14) 



(15) 



The functions ip(E, r) are routinely available in computations starting at near origin 
(r = 0) with 

[^{E,r)] n = r l 5 vl . 
If in asymptotic region, for r > S, these radial functions behave as 

i 



then 



[i>{E,r)] in =j v WEr)[C{E)] in 



[C{E)] in = W [iP(E, r)] in , —=h v (VEr) 



h v {^Er)[S{E)\ in 



(16) 
(17) 



E 



r=S 



[S(E)] in =W {ji'(VEr),[^(E,r)} in } r _ 
with the Wronskian definition 



(18) 
(19) 



W{f,g} = f§- r g-g§- r f. 



Comparing (|11|) and H17|) gives us the following relation for the i-matrix 

t(E) = S(E)C^ 1 (E) . 

The poles of GF in 11U|) determine energy dispersion curves in the form of bands 
These poles are usually found from zeros of the KKR-determinant 



(20) 
(21) 



detlr 1 ^) -B{E,k)\ = 0. 



(22) 
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Furthermore, we define the matrix of logarithmic derivatives 

D(E,r) = —mE,r)}ty(E,r)}- 1 . (23) 

or 

The D(E, r) matrix is not affected by changing of the regular solution normalization. 
Moreover, this matrix is symmetrical, i.e. 

D T (E,r) = D(E,r) (24) 

what can be proved by recalling that it satisfies equation of Ricatti type 

^-D(E,r)+D 2 {E,r) = v(r) - E (25) 
or 

with the symmetrical potential matrix, v T (r) = v(r). At the 'muffin-tin' boundary 
point (r = S), we can simply write 



D(E)=^(E,r)} 



[TpiEiST 1 . (26) 

r=S 

The log-derivative matrix D(E) is also directly related to the i-matrix by the following 
expression: 

[v P \E)]i>i =3v{^ES)[t-\E)] vl]l {jES) 



l =h v {VES) ] i{^S)8vi 
E 

with the definition of the pseudopotential amplitude matrix 



(27) 



where 



Vp{E) = D(E) — Dq(E) (28) 



D (E) = ±\j(VEr)] [j(VeS)}-\ 
ar r=s 



(29) 

\j(VEr)]vi =MVEt)S V i . 
3. Novel forms of the KKR-matrix in one dimension 

In the standard KKR-methodology finding zeros of the matrix [t~ 1 (E) — B(E, k)] in 
the GF formula I|1U|I is central topic when computing energy bands. This is done by 
the requirement that at least one of the eigenvalues of that matrix goes through zero. 
In fact that condition comes out from analyzing not the full GF in i|10|) but only from 
its second part. To circumvent this problem, we want the first term in Ea. H10|) to 
vanish at the boundary radial points. It can be done by re-normalizing the regular 
solution as 

£(£,r) =Z(E,r)Z- 1 (E,S) (30) 
and accordingly the irregular solution 

C(E, r) = J{E, r)Z T (E, S) - Z(E, r)j(VES). (31) 
By direct calculations we may check that the Wronskian in matrix form defined as 

W{aE,rte(E,r)} vl = jc 7 \E , r) ^{E \ r) - |-C T (£, r)£(E, r)| ^ (32) 
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is equal to 

W{((E, r),i(E, r)} = W{J(E, r),Z(E, r)} = 1. (33) 

Setting in (J30JI and lO r = S and recalling that t 7 '(E) = t(E) it follows, that 

£(E,S) = 1 ((E,S)=0. (34) 

Now, we are in position to convert GF expressed in using convention (J, Z) to 
that in terms of (C,£). First, we need the free GF matrix with radii on 'muffin-tin' 
points: 

[go{E, k)]in = -^=h l (VES)ji(VES)6 in 

^ E (35) 



+j l ,{VES)[B(E,k)] vl j l (VES) . 

Then applying the operator identity 

[A - B}- 1 = A' 1 + A' 1 [B- 1 - A' 1 ] _1 A' 1 (36) 

to the second term in (|f Ofl and replacing radial functions (J,Z) with functions (£, £) 
as introduced in 1)31) [I and (|31ll we get to the following expression of GF: 

< x'\G(E,k)\x >=^2Y l2 (x') <hr'\G(E,k)\hr>Y h (x) (37) 
hh 

with radial parts 



< h r>\G(E,k)\h r >= - £[C(£, » >)]w[£(£, r<" T 



Hi 



VI 



(38) 



It follows from (|38|l that GF with radial arguments at boundary points 

[g(E, k)] hh =< l 2 S\G(E, k)\h S > (39) 

can be found from algebraic equation of Dyson type 

g(E, k) = g (E, k) + g (E, k)v P (E)g(E, k) (40) 

with pseudopotential amplitude vp(E) as given in 1)28(1. The poles of g(E,k) are 
exactly the same as of the full GF < x'\G(E, k)\x > in QlOfl. suggesting alternative 
way to <|22|) of finding energy bands 

det\g- 1 (E,k)\=0. (41) 

Searching for zeros of eigenvalues of matrix g~ 1 (E,k) is much easier than that 
of t — 1 (£") — B(E,k) as the eigenvalue functions of the former are monotonically 
increasing with energy E. But some obscuring deficiency still persists as the slopes, 
at which these eigenvalues are crossing energy axis, are not fixed. This inconvenience 
can be avoided in the following way. Let us first integrate GF in the form l|38|l over 
the whole Wigner-Seitz cell 

a/2 a/2 

J dx < x\G(E,k)\x >= -J2 I dx Q(E,x)£i{E,x) 

-a/2 I -a/2 

(42) 

a/2 

+ fc)]l'J / dx^(E,xME,x). 

I'l -a/2 
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Integrating the first term in (|42H yields: 

a/2 

J dxb>{E,x)£t(E,x) = 

-a/2 

E hh fdr Y h (x)Y h (i)[C(S, r)] hV [£(E, r)] hl = (43) 
o 

f<MC T (£,r)£(£,r)]^ 
o 

From the Schrodinger equation we find that radial parts £ T (E, r) and the energy 
derivative -^^(E 7 r) = £(E,r) satisfy differential equations 

^C T (E,r)=C T (E,r)[v(r)-E] (44) 

and 

^i(E, r) = [v(r) - E]£(E, r) - £(E, r) (45) 

respectively. Multiplying l|4^|) to the right by £(E, r), then l|43|l to the left by C T (-Ej r ) 
and subtracting one from another, we get: 

s 

drC T (E,r)aE,r) = W {((E, r),i(E, r)}^ (46) 

This integral can be conventionally found if function ip{E, r) is computed. From 
the definition of £(E, r) in (|3U|l it follows: 

t(E,r)1>(E,S)=1>(E,r), (47) 

i(E,r)1>(E,S)+t(E,r)ij>(E,S)=i>(E,r) = (r - 0), (48) 

Inserting £(E,r) from (14811 into 14(i|) we get 

w{((E,r),i(E,r)\ = 

w\t(E,r),-t(E,r)i>(E,S)i>- l (E,S)} = (49) 



-i>(E,S)i/>-HE,S) 
Combining H46[) with the result (|49|l we have 



s 

dr( T (E, r)£{E, r) = —ip(E, S)^~ l (E, S) . (50) 



For the second integral in (|42ll we obtain 



T dx&(E,x)Z l (E ) x) = fdr[f{E,r)£(E,r)] vl = 
-o/ 2 (51) 

[N^E)^-\E, S)] T [N4E)^\E, S)] 
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with the definition of N^{E) as 
S 

dri> T (E,r)i>{E,r) = N^(E)N^(E) . (52) 

With the help of ifSUI) and (|5"T)l we may convert (|4"2|) to the form 

T < x\G(E,k)\x >= UHE^^-HE^S)],^ 

-a/2 I 

E{m,(E)^\E,S)]g(E,k)[N^E)^(E,S)] T } u = (53) 

i 

i 

The left hand side of equation (|53|l is equal to the trace of hermitian matrix P(E, k) 
defined as 

P(E,k) = l/2[Q(E,k) + Q\E,k)], (54) 

with 

Q(E,k) = ^(E,S)N- 1 (E)r 1 [^(E,S)N- 1 (E)} + 

(55) 

N^{E)^{E, S)g(E, k)[N^{E)^-\E, S)] T . 

The poles of the matrix P{E, k) give energy bands that can be found from zeros of 
the determinant 

detlP-^fc)! =0 . (56) 

From the analytical property of GF known in the literature as Herglotz property 10 
it is expected that eigenvalues of P~ 1 (E, fc)-matrix must increase monotonically with 
energy, i.e., 

■JL\[p-\E,k)]>0, (57) 

and with the slope equal to one if crossing energy axis. 

In the formula (|55l) the presence of the inverse matrix ip~ 1 (E,S) may cause 
numerical instability if matrix ip(E, S) is singular, what may happen in practice. 
To avoid such complications it is convenient to combine the first and second term in 
Ea. H54|) together with the following manipulation. First, we can compute integral in 
l\~>'2li using trick with energy derivative, like in l|44|) and (|45fl . Then with the known 
procedure we get 

S 

dr?p T (E,r)iP(E,r) = -W {ip(E,r),iP(E,r)} ^. (58) 
o 

In equation l|58|l the contribution at the origin point (r = 0) is set to be equal 
zero, what is justified with normalization of ?p(E,r) assumed in Ea. (|16[) . Recalling 
Wronskian definition, as set in l|20l) and log-derivatives in l|26(). we may proceed with 
getting 

rS , , . 

dr^ T {E,r)^{E 1 r)=^ T {E,S)[D^{E)-D^E)^{E 1 S) (59) 



(61) 
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From H59f) and (|52[1 it follows the equality 

[iP(E,S)N-\E)] T {d^E) - D^E)} [iP(E,S)N-\E)\ = 1, (60) 

useful in converting (|55|l to the desired form as follows 

Q(E,k) = [i>(E, S)Ntp(E)]~ 1 x 

{i>(E,S)N-\E) + g(E,kmE,S)N^(E)]-i T } 

From l|61|l with the help of equality H6U|) we may proceed to 
Q(E,k) = [ij(E,S)N^(E)}- 1 x 

(62) 

{[D 4 ,{E) -D^E)]' 1 +g(E,k)} [^{E,S)N^(E)]^ T . 

But recalling g(E, k) as given in (|40[1 we get 

[D^(E)-D^E)]- l +g(E,k) = 

\go\E,k) +Do(E) - D^E)}-\g^\E,k) + D Q (E) ~ D^(E)]x (63) 

[D^(E) - D^{E)]-\ 
that if inserted in (|62(l gives 
Q(E,k) = 

[i>(E, S)N~ 1 (E) + g (E, k)[D (E) - D^{E)]^{E, S)N~ 1 (E)y 1 x (64) 

y>(E,S)Nj\E) + g Q (E,k)[D (E) - D^E)]^(E,S)N^(E)} . 
With the Wronskian definition as in l|2(J|l last expression can be simplified further 

Q(E, k) = N^(E) {i>(E, S) + C(E, k)W T [^, j]}" 1 x 

(65) 

[i>{E, S) + C(E, k)W T [iPJ] } N- 1 (E) 
with the matrix C(E, k) formed from KKR structure functions B(E, k): 

[C(E > k)] n = -^=h l ,(VES)5 l , l + j l ,(VES)[B(E,k)}n (66) 
V E 



These two last equations give the most convenient way to construct the novel 
KKR-matrix P(E, k) with desired properties of the eigenvalues. 

4. Generalised Lloyd formula 

The Lloyd formula is often used to calculate the number of states for a 'muffin-tin' 
model potential. We derive it here in the form similar to that used in the paper by 
Kaprzyk and Bansil 

It is well known that the total number of states N(E) is done as 

N(E) = --Im[ dE V / dx < x\G{E,k)\x > . (67) 



keBZ 
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We show that the space integral in i|67|) can be expressed as a perfect energy derivative 
and the energy integration can be done explicitly. We start with (|53|l putting it in the 
form 

a/2 

/ dx < x\G{E, k)\x >= Tr{ip-\E, S)ip(E, S)}+ 
-o/ 2 (68) 

Tr{g(E, k)[N^{E)^-\E, S)} T [N^{E)^-\E, S)}} , 
then with the help of l|59|l we find (at r = S) 

N^{E)N il {E) = -V{E, S)b^,{E)i){E, S). (69) 
With the last result, can be written as 



a/2 

J dx < x\G(E,k)\x >-- 

-a/2 



(70) 



Tr 



{^{E, S)4>(E, S)}-Tt {g(E, k)D^(E)} . 
is of proving the equality 

g(E,k)=g(E,k)D^(E)g(E > k) (71) 



The last step consists of proving the equality 
d_ 
dE 

which if inserted into (|70|l gives the result 

a/2 d 

dx < x\G(E,k)\x >= -— Trln[ ff (S,ifc)i/'" 1 ( J B,5)]. (72) 

-a/2 OE 

The algebraic equality 171|) is a consequence of GF properties stating that 

-^G(E) = —G{E)G(E), (73) 

If applied to our Bloch-Fourier GF with radial arguments on 'muffin-tin' boundary 
points we get 

<x 2 S\-§ E G{E,k)\x 1 S> = 

a/2 (74) 

- j dx < x 2 S\G(E, k)\x >< x\G(E, k)\x! S > 

-a/2 

and for the radial part defined in i|37|) 

<l 2 S\J^G{E,k)\hS> = 



s 

Jdr< l 2 S\G(E, k)\l r >< I r\G(E, k)\h S > 
o 



(75) 



Inserting proper radial arguments into (|75|l we get according to i|36|) and (|34[l that 

< l 2 S\G(E, k)\l r >= [g(E, k)( T (E, r)] hl , (76) 

< I S\G(E, k)\h r >= [£(E, r)g(E, k)} lh , (77) 
which allow us to rewrite l|75|) in the form 



-^[g(E,k)] = -g(E,k) 



J drf(E,r)Z(E, 



g(E,k). (78) 
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Computing integral in (|75|l can be done with the help of and leads to the result 

S 

J drf{E, r)i{E, r) = -D 6 (E) = -D^E). (79) 
o 

Inserting this relation into H78[) proves validity of (|71|l . The total number of states 
seen in (I67f) can now be found from 

N(E) = -lmln{det\Tp(E,S)\) - - V Imlndet \g(E, k)\. (80) 

k£BZ 

In passing from Ij72(l to (|80|l we used algebraic equality 

Trln[A] = ln[det|A|]. (81) 

In l|80l) the first term accounts for the number of nodes on regular solution from origin 
up to radius S, and the second term comes out from the Bloch states at each of 
/c-points in BZ. 

5. Results 

In order to illustrate numerically how the novel approach to KKR method works 
compared to standard one, we have performed calculations for the ID Mathieu 
potential of the form 

v{x) = -U cos(2irx/a), (82) 

and used lattice constant a — 3.0 a.u. and Uq — 5 Ry (see Fig. In Fig. we show 
eigenvalue functions as calculated using both standard and novel method (Eqs. (I22|l 
and ((23)) for the potential l|82|l . One can see that zeros of the eigenvalue functions are 
exactly the same. In the novel method case we can find them numerically very easily 
but in the standard KKR approach it is much more difficult. Moreover, to illustrate 
how the eigenvalue curves behave when the number of atoms in the unit cell grows, we 
performed the calculations in two cases. In the first case we put only one atom in the 
unit cell with the 'muffin-tin'' radius S = 1.5 a.u. and the symmetric potential given 
by H82|) . In the second case we divided the unit cell into four 'muffin-tin' segments 
of the same length and touching each other, but without changing the total potential 
in the unit cell. So, in this case the potential in each segment is not symmetrical, 
but the physical situation is exactly the same. What changed is the size of resulting 
matrices and hence number of the eigenvalue functions. These functions in both cases 
are shown in Fig.[3|for the wave vector k = 0.67r/a. Note that in the second case the 
lines are not derived from interpolation but simply connect calculated points, which 
are not shown in Fig. [3] 

Now all the benefits of the novel method can be underlined. First, we observe that 
the eigenvalue functions do not cross each other and grow monotonically. Moreover, 
they form almost straight lines with the slope of unity when passing through zero. 
Secondly, we can see that with the increasing number of atoms in the unit cell the 
eigenvalue curves become straight for the wider range of energies. So, in the limit of 
infinite number of atoms in the unit cell, the method seems to become linear. This 
makes the method not very sensitive to the number of energy points used in the 
interpolation procedure, when finding zeros of eigenvalue curves. This last advantage 
opens possibilities of calculation of very complex systems, even with hundreds atoms 
in the unit cell in reasonable time without loosing accuracy. 
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-6 1 1 ■ 1 1 1 1 1 

-1.5 -1 -0.5 0.5 1 1.5 

x [a.u.] 

Figure 1. Mathieu potential used in calculations. The black dots denote four 
atoms in the unit cell case. In the case of one atom, it is placed in the middle of 
the unit cell. 



6. Summary 

In this paper we have reformulated one-dimensional full potential KKR formalism. 
Then we have derived novel form of the KKR secular matrix. To exclude first term in 
the full Green function we have normalized solutions of the Schrodinger equation in 
such a way that the 'irregular' solution disappears at boundary points of 'muffin-tin' 
segments. Then by integrating the full Green function over the Wigner-Seitz cell we 
have derived the final result, that is the expressions for the secular matrix P(E, k) 
as well as for the generalised Lloyd formula. The eigenvalues of the inverse matrix 
P^ 1 (E,k) increase monotonically with energy and almost linearly for every /c-point. 
In the case of increasing number of atoms in the unit cell the eigenvalue functions 
become more linear. 

Finally, we performed numerical calculations for the case of the Mathieu potential. 
The results show that the zeros of the eigenvalues can be easily found, even for very 
complex systems, without any lost of accuracy. This formalism can be extended to 
the higher-dimensional systems with the minor changes. 

Appendix 

In this appendix we briefly generalise the formalism introduced in sections 2, 3 and 4 
for many atoms in the unit cell. In this case one has to solve the Schrodinger equation 
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Figure 2. Comparison of the standard KKR method (dots) and the novel 
one (solid lines) as applied to the Mathieu potential 1821 . The wave vector is 
k = 0.67r/a. 



with the potential 

oo p 

V(x) = 2J ^^(l-M-fla), (A.l) 
n— — oo a— 1 

where a Q is the position of a-th atom in the unit cell with p equal to the number 

of atoms in the unit cell. At each 'muffin-tin' atomic segment we solve Schrodinger 
equation starting at origin from 

[i> a {E,r)] vl = r l 5 vl r^O. (A.2) 

The resulting ip matrix has the form 

bP(E, r)] a , v , a i = [tp a (E, r)] in 5 aa > (A.3) 

with the size 2p x 2p. Similarly we construct all other matrices introduced in the 
previous sections. GF in (£,£) representation now has a form: 

< x' + a a >\G(E, k)\x + a a >= 



Y,Qi.E,x > )^{E 1 x < )8 ac 



■ £ (E, aO [So \E,k)-v P (E)]-}, (E, x) 
1,1' 



(AA) 
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Figure 3. The eigenvalue functions for the Mathieu potential 1821 and the wave 
vector k = O.Gir/a. The case with one and four atoms in the unit cell is represented 
by solid and dashed lines, respectively. 



The structure constants in this generalised case can be given explicitly with the matrix 
elements [B(E,k)] a >i',aV 



[B(E,k)] a o,a'o = [B{E,k)} 

exp(iy/Ea) cos^fEa aa ' — cos(fca — \f~Ea 
iyE[cos ka — cos s/Ea] 



(A.5) 



H 7= exp(«V£;|cw|)(l - 6 aa >), 

IV E 

[B(E,k)U, a , = [B(E,k)}* a , 0tal = 

exp^V^Ea) sin \f~Ea aa i + sin(fca — y/Eaaa 1 ) 
i\fE[co'& ka — cos v^Ea] 

H — j=exp(i- v / B|a aQ /|)sgn(-a Q!a ,')(l - S aa >), 
\E 



(A.6) 



with a aa i — a a — a a i. Finally we define the matrix P(E, k) by the relation: 

Sc. 

^2 / dx<x + a a \G(E,k)\x + a a >=Tr[P(E,k)] (A.7) 
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with S a being muffin-tin radius of the a-th atom. Following the steps described in 
previous sections we find, that the matrix P(E, k) can be found using Eqs. (|54H 
and l|55|) by adding atomic indexes to the matrix elements, as it is done for structure 
constants matrix B(E, k). The band structure is still determined by the relation 
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